Genomic Diversity of Bacteriophages Infecting the Genus Acinetobacter.

The number of sequenced Acinetobacter phage genomes in the International Nucleotide Sequence Database Collaboration has increased significantly in recent years, from 37 in 2017 to a total of 139 as of January 2021 with genome sizes ranging from 31 to 378 kb. Here, we explored the genetic diversity of the Acinetobacter phages using comparative genomics approaches that included assessment of nucleotide similarity, shared gene content, single gene phylogeny, and the network-based classification tool vConTACT2. Phages infecting Acinetobacter sp. are genetically diverse and can be grouped into 8 clusters (subfamilies) and 46 sub-clusters (genera), of which 8 represent genomic singletons (additional genera). We propose the creation of five new subfamilies and suggest a reorganisation of the genus Obolenskvirus. These results provide an updated view of the viruses infecting Acinetobacter species, providing insights into their diversity.

Genome Analysis and Therapeutic Evaluation of a Novel Lytic Bacteriophage of Salmonella Typhimurium: Suggestive of a New Genus in the Subfamily Vequintavirinae.

Salmonella Typhimurium, a foodborne pathogen, is a major concern for food safety. Its MDR serovars of animal origin pose a serious threat to the human population. Phage therapy can be an alternative for the treatment of such MDR Salmonella serovars. In this study, we report on detailed genome analyses of a novel Salmonella phage (Salmonella-Phage-SSBI34) and evaluate its therapeutic potential. The phage was evaluated for latent time, burst size, host range, and bacterial growth reduction in liquid cultures. The phage stability was examined at various pH levels and temperatures. The genome analysis (141.095 Kb) indicated that its nucleotide sequence is novel, as it exhibited only 1-7% DNA coverage. The phage genome features 44% GC content, and 234 putative open reading frames were predicted. The genome was predicted to encode for 28 structural proteins and 40 enzymes related to nucleotide metabolism, DNA modification, and protein synthesis. Further, the genome features 11 tRNA genes for 10 different amino acids, indicating alternate codon usage, and hosts a unique hydrolase for bacterial lysis. This study provides new insights into the subfamily Vequintavirinae, of which SSBI34 may represent a new genus.

A Polyvalent Broad-Spectrum Escherichia Phage Tequatrovirus EP01 Capable of Controlling Salmonella and Escherichia coli Contamination in Foods.

Salmonella and Escherichia coli (E. coli) food contamination could lead to serious foodborne diseases. The gradual increase in the incidence of foodborne disease invokes new and efficient methods to limit food pathogenic microorganism contamination. In this study, a polyvalent broad-spectrum Escherichia phage named Tequatrovirus EP01 was isolated from pig farm sewage. It could lyse both Salmonella Enteritidis (S. Enteritidis) and E. coli and exhibited broad host range. EP01 possessed a short latent period (10 min), a large burst size (80 PFU/cell), and moderate pH stability (4-10) and appropriate thermal tolerance (30-80 °C). Electron microscopy and genome sequence revealed that EP01 belonged to T4-like viruses genus, Myoviridae family. EP01 harbored 12 CDSs associated with receptor-binding proteins and lacked virulence genes and drug resistance genes. We tested the inhibitory effect of EP01 on S. Enteritidis, E. coli O157:H7, E. coli O114:K90 (B90), and E. coli O142:K86 (B) in liquid broth medium (LB). EP01 could significantly reduce the counts of all tested strains compared with phage-free groups. We further examined the effectiveness of EP01 in controlling bacterial contamination in two kinds of foods (meat and milk) contaminated with S. Enteritidis, E. coli O157:H7, E. coli O114:K90 (B90), and E. coli O142:K86 (B), respectively. EP01 significantly reduced the viable counts of all the tested bacteria (2.18-6.55 log10 CFU/sample, p < 0.05). A significant reduction of 6.55 log10 CFU/cm2 (p < 0.001) in bacterial counts on the surface of meat was observed with EP01 treatment. Addition of EP01 at MOI of 1 decreased the counts of bacteria by 4.3 log10 CFU/mL (p < 0.001) in milk. Generally, the inhibitory effect exhibited more stable at 4 °C than that at 28 °C, whereas the opposite results were observed in milk. The antibacterial effects were better at MOI of 1 than that at MOI of 0.001. These results suggests that phage EP01-based method is a promising strategy of controlling Salmonella and Escherichia coli pathogens to limit microbial food contamination.

Isolation, characterization, molecular analysis and application of bacteriophage DW-EC to control Enterotoxigenic Escherichia coli on various foods.

Among food preservation methods, bacteriophage treatment can be a viable alternative method to overcome the drawbacks of traditional approaches. Bacteriophages are naturally occurring viruses that are highly specific to their hosts and have the capability to lyse bacterial cells, making them useful as biopreservation agents. This study aims to characterize and determine the application of bacteriophage isolated from Indonesian traditional Ready-to-Eat (RTE) food to control Enterotoxigenic Escherichia coli (ETEC) population in various foods. Phage DW-EC isolated from Indonesian traditional RTE food called dawet with ETEC as its host showed a positive result by the formation of plaques (clear zone) in the bacterial host lawn. Transmission electron microscopy (TEM) results also showed that DW-EC can be suspected to belong to the Myoviridae family. Molecular characterization and bioinformatic analysis showed that DW-EC exhibited characteristics as promising biocontrol agents in food samples. Genes related to the lytic cycle, such as lysozyme and tail fiber assembly protein, were annotated. There were also no signs of lysogenic genes among the annotation results. The resulting PHACTS data also indicated that DW-EC was leaning toward being exclusively lytic. DW-EC significantly reduced the ETEC population (P ≤ 0.05) in various food samples after two different incubation times (1 day and 6 days) in chicken meat (80.93%; 87.29%), fish meat (63.78%; 87.89%), cucumber (61.42%; 71.88%), tomato (56.24%; 74.51%), and lettuce (46.88%; 43.38%).

Characterization and Genome Study of Novel Lytic Bacteriophages against Prevailing Saprophytic Bacterial Microflora of Minimally Processed Plant-Based Food Products.

The food industry is still searching for novel solutions to effectively ensure the microbiological safety of food, especially fresh and minimally processed food products. Nowadays, the use of bacteriophages as potential biological control agents in microbiological food safety and preservation is a promising strategy. The aim of the study was the isolation and comprehensive characterization of novel bacteriophages with lytic activity against saprophytic bacterial microflora of minimally processed plant-based food products, such as mixed leaf salads. From 43 phages isolated from municipal sewage, four phages, namely Enterobacter phage KKP 3263, Citrobacter phage KKP 3664, Enterobacter phage KKP 3262, and Serratia phage KKP 3264 have lytic activity against Enterobacter ludwigii KKP 3083, Citrobacter freundii KKP 3655, Enterobacter cloacae KKP 3082, and Serratia fonticola KKP 3084 bacterial strains, respectively. Transmission electron microscopy (TEM) and whole-genome sequencing (WGS) identified Enterobacter phage KKP 3263 as an Autographiviridae, and Citrobacter phage KKP 3664, Enterobacter phage KKP 3262, and Serratia phage KKP 3264 as members of the Myoviridae family. Genome sequencing revealed that these phages have linear double-stranded DNA (dsDNA) with sizes of 39,418 bp (KKP 3263), 61,608 bp (KKP 3664), 84,075 bp (KKP 3262), and 148,182 bp (KKP 3264). No antibiotic resistance genes, virulence factors, integrase, recombinase, or repressors, which are the main markers of lysogenic viruses, were annotated in phage genomes. Serratia phage KKP 3264 showed the greatest growth inhibition of Serratia fonticola KKP 3084 strain. The use of MOI 1.0 caused an almost 5-fold decrease in the value of the specific growth rate coefficient. The phages retained their lytic activity in a wide range of temperatures (from -20 °C to 50 °C) and active acidity values (pH from 4 to 11). All phages retained at least 70% of lytic activity at 60 °C. At 80 °C, no lytic activity against tested bacterial strains was observed. Serratia phage KKP 3264 was the most resistant to chemical factors, by maintaining high lytic activity across a broader range of pH from 3 to 11. The results indicated that these phages could be a potential biological control agent against saprophytic bacterial microflora of minimally processed plant-based food products.

Genome and Ecology of a Novel Alteromonas Podovirus, ZP6, Representing a New Viral Genus, Mareflavirus.

Alteromonas is a ubiquitous, abundant, copiotrophic and phytoplankton-associated marine member of the Gammaproteobacteria with a range extending from tropical waters to polar regions and including hadal zones. Here, we describe a novel Alteromonas phage, ZP6, that was isolated from surface coastal waters of Qingdao, China. ZP6 contains a linear, double-stranded, 38,080-bp DNA molecule with 50.1% G+C content and 47 putative open reading frames (ORFs). Three auxiliary metabolic genes were identified, encoding metal-dependent phosphohydrolase, diaminopurine synthetase, and nucleotide pyrophosphohydrolase. The first two ORFs facilitate the replacement of adenine (A) by diaminopurine (Z) in phage genomes and help phages to evade attack from host restriction enzymes. The nucleotide pyrophosphohydrolase enables the host cells to stop programmed cell death and improves the survival rate of the host in a nutrient-depleted environment. Phylogenetic analysis based on the amino acid sequences of whole genomes and comparative genomic analysis revealed that ZP6 is most closely related to Enhodamvirus but with low similarity (shared genes, <30%, and average nucleotide sequence identity, <65%); it is distinct from other bacteriophages. Together, these results suggest that ZP6 could represent a novel viral genus, here named Mareflavirus. Combining its ability to infect Alteromonas, its harboring of a diaminopurine genome-biosynthetic system, and its representativeness of an understudied viral group, ZP6 could be an important and novel model system for marine virus research. IMPORTANCEAlteromonas is an important symbiotic bacterium of phytoplankton, but research on its bacteriophages is still at an elementary level. Our isolation and genome characterization of a novel Alteromonas podovirus, ZP6, identified a new viral genus of podovirus, namely, Mareflavirus. The ZP6 genome, with a diaminopurine genome-biosynthetic system, is different from those of other isolated Alteromonas phages and will bring new impetus to the development of virus classification and provide important insights into novel viral sequences from metagenomic data sets.

Isolation and characterization of a lytic bacteriophage against Pseudomonas aeruginosa.

In recent years, the use of bacteriophages (or 'phages') against multidrug-resistant (MDR) bacteria including Pseudomonas aeruginosa has drawn considerable attention, globally. In this work, we report the isolation and detailed characterization of a highly lytic Pseudomonasphage DRL-P1 isolated from wastewater. Under TEM, DRL-P1 appeared as a member of the phage family Myoviridae. DRL-P1 featured rapid adsorption (~ 5 min), short-latency (~ 30 min), and large burst size (~ 100 PFU per infected cell). DRL-P1 can withstand a wide temperature range (4 °C to 40 °C) and pH (5.0 to 10.0) conditions. The 66,243 bp DRL-P1 genome (MN564818) encodes at least 93 ORFs, of which 36 were functionally annotated based on homology with similar phage proteins available in the databases. Comparative analyses of related genomes suggest an independent evolutionary history and discrete taxonomic position of DRL-P1 within genus Pbunavirus. No toxin or antibiotic resistance genes was identified. DRL-P1 is tolerant to lyophilization and encapsulation techniques and retained lytic activity even after 18 months of storage. We also demonstrated decontaminating potentials of DRL-P1 in vitro, on an artificially contaminated cover-slip model. To the best of our knowledge, this is the first Pbunavirus to be reported from India. Our study suggests DRL-P1 as a potential candidate for various applications.

Novel Viruses That Lyse Plant and Human Strains of Kosakonia cowanii.

Kosakonia cowanii (syn. Enterobacter cowanii) is a highly competitive bacterium that lives with plant, insect, fish, bird, and human organisms. It is pathogenic on some plants and an opportunistic pathogen of human. Nine novel viruses that lyse plant pathogenic strains and/or human strains of K. cowanii were isolated, sequenced, and characterized. Kc166A is a novel kayfunavirus, Kc261 is a novel bonnellvirus, and Kc318 is a new cronosvirus (all Autographiviridae). Kc237 is a new sortsnevirus, but Kc166B and Kc283 are members of new genera within Podoviridae. Kc304 is a new winklervirus, and Kc263 and Kc305 are new myoviruses. The viruses differ in host specificity, plaque phenotype, and lysis kinetics. Some of them should be suitable also as pathogen control agents.

Characterization and genome sequence of the genetically unique Escherichia bacteriophage vB_EcoM_IME392.

In this study, a novel Escherichia coli-specific bacteriophage, vB_EcoM_IME392, was isolated from chicken farm sewage in Qingdao, China. The genome of IME392 was found by next-generation sequencing to be 116,460 base pairs in length with a G+C content of 45.4% (GenBank accession number MH719082). BLASTn results revealed that only 2% of the genome sequence of IME392 shows sequence similarity to known phage sequences in the GenBank database, which indicates that IME392 is a novel bacteriophage. Transmission electron microscopy showed that IME392 belongs to the family Myoviridae. The host range, the multiplicity of infection, and a one-step growth curve were also determined.

Incidence and characterisation of Bacillus cereus bacteriophages from Thua Nao, a Thai fermented soybean product.

Bacillus cereus is considered to be an important food poisoning agent causing diarrhea and vomiting. In this study, the occurrence of B. cereus bacteriophages in Thai fermented soybean products (Thua Nao) was studied using five B. cereus sensu lato indicator strains (four B. cereus strains and one B. thuringiensis strain). In a total of 26 Thua Nao samples, there were only two bacteriophages namely BaceFT01 and BaceCM02 exhibiting lytic activity against B. cereus. Morphological analysis revealed that these two bacteriophages belonged to the Myoviridae. Both phages were specific to B. cereus and not able to lyse other tested bacteria including B. licheniformis and B. subtilis. The two phages were able to survive in a pH range between 5 and 12. However, both phages were inactive either by treatment of 50°C for 2 h or exposure of UV for 2 h. It should be noted that both phages were chloroform-insensitive, however. This is the first report describing the presence of bacteriophages in Thua Nao products. The characterization of these two phages is expected to be useful in the food industry for an alternative strategy including the potential use of the phages as a biocontrol candidate against foodborne pathogenic bacteria.

Complete genome analysis of the newly isolated Shigella sonnei phage vB_SsoM_Z31.

This work describes the characterization and genome annotation of the newly isolated lytic phage vB_SsoM_Z31 (referred to as Z31), isolated from wastewater samples collected in Dalian, China. Transmission electron microscopy revealed that phage Z31 belongs to the family Myoviridae, order Caudovirales. This phage specifically infects Shigella sonnei, Shigella dysenteriae, and Escherichia coli. The genome of the phage Z31 is an 89,355-bp-long dsDNA molecule with a G+C content of 38.87%. It was predicted to contain 133 ORFs and encode 24 tRNAs. No homologs of virulence factor genes or antimicrobial resistance genes were found in this phage. Based on the results of nucleotide sequence alignment and phylogenetic analysis, phage Z31 was assigned to the genus Felixounavirus, subfamily Ounavirinae.

Characterization of a lytic EBP bacteriophage with large size genome against Enterobacter cloacae.

Enterobacter cloacae (E. cloacae) is an emerging nosocomial pathogen that had acquired antibiotic resistance against multiple classes of antibiotics. The current study was aimed to isolate and characterize lytic bacteriophage against E. cloacae. The bacteriophage EBP was isolated from a sewage water sample using E. cloacae as a host strain by double-layer agar technique. EBP was found stabile at a wide range of temperatures (25, 37, 60, and 80°C) and pH (5, 6, 7, 8, and 9) with antibacterial activity up to 24 h of infection. The latent period of EBP was 20 min with a burst size of 252 phages per cell. It showed a narrow host range and infected 12/21 (57%) isolates of E. cloacae tested. It has helical symmetry with a head size of 105 and 120 nm long tail with contractile sheath. The EBP has 179.1 kb long double-stranded DNA genome with 44.8% GC content. Majority of identified ORFs (187/281) were encoding putative proteins with unknown function. Necessary replication enzymes, structural proteins, and lytic enzymes were detected in the genome of EBP. Phylogenetic analysis revealed that EBP closely resembles with Coronobacter phage vB_CsaM_IeN, vB_CsaM_IeE, vB_CsaM_IeB, and Citrobacter phage Margaery. Based on electron microscopy and molecular characterization, EBP was classified as a Myoviridae phage.

First Description of a Temperate Bacteriophage (vB_FhiM_KIRK) of Francisella hispaniensis Strain 3523.

Here we present the characterization of a Francisella bacteriophage (vB_FhiM_KIRK) including the morphology, the genome sequence and the induction of the prophage. The prophage sequence (FhaGI-1) has previously been identified in F. hispaniensis strain 3523. UV radiation induced the prophage to assemble phage particles consisting of an icosahedral head (~52 nm in diameter), a tail of up to 97 nm in length and a mean width of 9 nm. The double stranded genome of vB_FhiM_KIRK contains 51 open reading frames and is 34,259 bp in length. The genotypic and phylogenetic analysis indicated that this phage seems to belong to the Myoviridae family of bacteriophages. Under the conditions tested here, host cell (Francisella hispaniensis 3523) lysis activity of KIRK was very low, and the phage particles seem to be defective for infecting new bacterial cells. Nevertheless, recombinant KIRK DNA was able to integrate site-specifically into the genome of different Francisella species after DNA transformation.

Characteristics of Bacteriophage Isolates and Expression of Shiga Toxin Genes Transferred to Non Shiga Toxin-Producing E. coli by Transduction.

A risk analysis of Shiga toxin (Stx)-encoding bacteriophage was carried out by confirming the transduction phage to non-Stx-producing Escherichia coli (STEC) and subsequent expression of the Shiga toxin genes. The virulence factor stx1 was identified in five phages, and both stx1 and stx2 were found in four phages from a total of 19 phage isolates with seven non-O157 STEC strains. The four phages, designated as φNOEC41, φNOEC46, φNOEC47, and φNOEC49, belonged morphologically to the Myoviridae family. The stabilities of these phages to temperature, pH, ethanol, and NaClO were high with some variabilities among the phages. The infection of five non-STEC strains by nine Stx-encoding phages occurred at a rate of approximately 40%. Non-STEC strains were transduced by Stx-encoding phage to become lysogenic strains, and seven convertant strains had stx1 and/or stx2 genes. Only the stx1 gene was transferred to the receptor strains without any deletion. Gene expression of a convertant having both stx1 and stx2 genes was confirmed to be up to 32 times higher for Stx1 in 6% NaCl osmotic media and twice for Stx2 in 4% NaCl media, compared with expression in low-salt environments. Therefore, a new risk might arise from the transfer of pathogenic genes from Stx-encoding phages to otherwise harmless hosts. Without adequate sterilization of food exposed to various environments, there is a possibility that the toxicity of the phages might increase.

Genomic analysis of bacteriophage Xoo-sp13 infecting Xanthomonas oryzae pv. oryzae.

Xanthomonas oryzae pv. oryzae is a bacterial pathogen that gives rise to diseases in rice all over the world. A bacteriophage infecting this bacterium was isolated from rice fields in China. Here, we report the complete genome sequence of this phage, which has a linear dsDNA genome of 309,023 bp and a G + C content of 42.43%. It contains 401 open reading frames and encodes 28 tRNAs. It belongs to the family Myoviridae and has a broad host range, making it a possible candidate for phage therapy.

Isolation, characterization and genomic analysis of vB-AhyM-AP1, a lytic bacteriophage infecting Aeromonas hydrophila.

AIMS: Aeromonas hydrophila is a zoonotic pathogen displaying resistance to multiple antibiotics. Here, we aim to develop a candidate biocontrol agent against A. hydrophila. METHODS AND RESULTS: In this study, we isolated and characterized the phage vB-AhyM-AP1 from sewage. It showed lytic activity against A. hydrophila strains. One-step growth curve revealed that the latent period lasted for 40 min. The burst size of one lytic cycle was 1413 PFU per infected cell. Temperature stability studies showed that the phage vB-AhyM-AP1 was active over temperatures ranging from 4 to 45°C for 1 h. pH stability studies indicated that the phage remained active within a pH range of 5-10 after 24 h of incubation. Stability tests in salt solutions showed that the phage was stable at salinities ranging from 0·1 to 2%. The phage also showed stabilities in organic solvents when incubated for 10 min. The Illumina Hiseq sequencing of its genome indicated that the phage vB-AhyM-AP1was a jumbo phage with a genome size of 2, 54 490 bp and GC content of 40·3%. The phylogenetic analysis of the terminase large subunit and major capsid protein indicated that the phage closely clustered with other Tevenvirinae phages. The genome encoded 455 ORFs and 22 tRNAs. The phage resulted in a reduction of 0·8 log units of viable A. hydrophila cells in biofilms grown on PVC coupons maintained in a low nutrient medium for 10 days. CONCLUSIONS: The phage showed lytic activity against planktonic and biofilm cells of A. hydrophila. Genome-based prediction showed it to be a strictly lytic phage without any virulence or antibiotic resistance genes indicating safety for environmental and clinical applications. SIGNIFICANCE AND IMPACT OF THE STUDY: The multidrug-resistant strains of A. hydrophila pose a significant health risk to both cultured fishes and consumers leaving few options for treatment. Phage vB-AhyM-AP1 may be used as a candidate biocontrol agent against A. hydrophila strains.

Genome analysis of Salmonella phage vB_SalM_8-19 (genus Rosemountvirus).

This work describes the characterization and genome annotation of Salmonella phage vB_SalM_8-19 (referred to as 8-19) isolated from sewage samples collected in a pig farm in Jilin, China. This phage was capable of infecting 60% Salmonella strains in our lab stock. The genome of phage 8-19 is composed of linear double-stranded DNA that is 52,648 bp in length with a G + C content of 46.02%; containing 74 ORFs and no tRNA genes. In October 2019, phylogenetic analyses indicated that phage 8-19 might belong to a novel cluster among the other similar phages which have not been specifically classified within some new genus in family Myoviridae. Recently, the International Committee on Taxonomy of Viruses (ICTV) defined phage 8-19 and its related phages as genus Rosemountvirus, family Myoviridae. This new genus, known as Rosemountvirus, is rarely reported in the literature.

T4-like myovirus community shaped by dispersal and deterministic processes in the South China Sea.

As the most abundant and genetically diverse biological entities, viruses significantly influence ecological, biogeographical and evolutionary processes in the ocean. However, the biogeography of marine viruses and the drivers shaping viral community are unclear. Here, the biogeographic patterns of T4-like viruses and the relative impacts of deterministic (environmental selection) and dispersal (spatial distance) processes were investigated in the northern South China Sea. The dominant viral operational taxonomic units were affiliated with previously defined Marine, Estuary, Lake and Paddy Groups. A clear viral biogeographic pattern was observed along the environmental gradient from the estuary to open sea. Marine Groups I and IV had a wide geographical distribution, whereas Marine Groups II, III and V were abundant in lower-salinity continental or eutrophic environments. A significant distance-decay pattern was noted for the T4-like viral community, especially for those infecting cyanobacteria. Both deterministic and dispersal processes influenced viral community assembly, although environmental selection (e.g. temperature, salinity, bacterial abundance and community, etc.) had a greater impact than spatial distance. Network analysis confirmed the strong association between viral and bacterial community composition, and suggested a diverse ecological relationship (e.g. lysis, co-infection or mutualistic) between and within viruses and their potential bacterial hosts.

Characterization of Novel Erwinia amylovora Jumbo Bacteriophages from Eneladusvirus Genus.

Jumbo phages, which have a genome size of more than 200 kb, have recently been reported for the first time. However, limited information is available regarding their characteristics because few jumbo phages have been isolated. Therefore, in this study, we aimed to isolate and characterize other jumbo phages. We performed comparative genomic analysis of three Erwinia phages (pEa_SNUABM_12, pEa_SNUABM_47, and pEa_SNUABM_50), each of which had a genome size of approximately 360 kb (32.5% GC content). These phages were predicted to harbor 546, 540, and 540 open reading frames with 32, 34, and 35 tRNAs, respectively. Almost all of the genes in these phages could not be functionally annotated but showed high sequence similarity with genes encoded in Serratia phage BF, a member of Eneladusvirus. The detailed comparative and phylogenetic analyses presented in this study contribute to our understanding of the diversity and evolution of Erwinia phage and the genus Eneladusvirus.

Susceptibility of Pseudomonas aeruginosa veterinary isolates to Pbunavirus PB1-like phages.

In recent years, antimicrobial-resistant Pseudomonas aeruginosa strains have increased in the veterinary field. Therefore, phage therapy has received significant attention as an approach for overcoming antimicrobial resistance. In this context, we isolated and characterized four Pseudomonas bacteriophages. Phylogenetic analysis showed that the isolated phages are novel Myoviridae Pbunavirus PB1-like phages with ØR12 belonging to a different clade compared with the other three. These phages had distinct lytic activity against 22 P. aeruginosa veterinary isolates. The phage cocktail composed from the PB1-like phages clearly inhibited the occurrence of the phage-resistant variant, suggesting that these phages could be useful in phage therapy.